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ABSTRACT 


Supersonic flow past oscillating cylindrical shells and 
oscillating flat~plate cascades is analyzed using linearized 
characteristics methods. Pressure distributions and gener- 
elazed aerodynamic forces are computed for arbitrary radius 
to length ratios, axial and circumferential wave numbers, 
Mach numbers, and reduced frequencies. Excellent agreement 
is obtained with a previous solution of the complete un- 
steady linearized potential equation by Dowell and Widnall 
using Laplace transform technigues. Also, pressure dis- 
tributions are computed for flat-plate cascades with sub- 
sonic leading edge locus for arbitrary solidity, stagger 
angle, frequency, and interblade phase angle. For compar- 
ison, a two-blade solution is developed using the method 
Of singularities. Good agreement is obtained for this 


Special case. 
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I. INTRODUCTION 


miecrsonte L£Liutter of skin panels first occurred on the 
German V-2 missile and has become of increasing engineering 
importance in recent years with the development of super- 
sonic flight vehicles. Panel flutter, for example, was 
found to occur on the X-15 Rocket Aircraft and Saturn Space 
Launch Vehicle. An integral part of the flutter analysis of 
such panels is the determination of the oscillatory pressure 
distribution. A number of approaches have been formulated 
Pemmen resulted mostly in approximation theories. A detailed 
review of the state of the art can be found in Reference 6. 
Only one “exact" approach was worked out by Dowell and Widnail 
using the complete linearized unsteady potential equation 
(Ref. 2). However, this approach requires rather complicated 
and time-consuming complex integration techniques. 

In the opesent work a method of characteristics procedure 
is developed to determine supersonic flow past vibrating 
panels and mein cal shells resulting in a rapid computation 
of the pressure distribution and generalized aerodynamic 
Forces. 

A second flutter problem has also become of considerable 
interest very recently. Due to the continuing trend toward 
higher tip speeds in modern jet engine fans and compressors, 
the possibility of supersonic blade flutter has become of 


considerable concern. A recent review of this problem area 


Je 





has been given in "A Workshop on Aeroelasticity in Turbo- 
machines" (Ref. 13). 

In the present study the method of characteristics is 
applied to compute supersonic flow past a finite flat plate 
cascade with subsonic leading edge locus. In the past, cas- 
cade studies have proved to be a useful tool for the under- 
standing of the complicated three-dimensional, incompressible 
or subsonic flow through an actual compressor stage. Al- 
though it is doubtful that the much more complex mixed 
transonic-~supersonic flow can be modeled by a two-dimensional 
cascade idealization, the analysis of an infinitely thin flat 
plate cascade may point the way toward a fuller understanding. 
The method of characteristics again proved to be a convenient 
tool for the analysis of supersonic fiow past oscillating 


cascades. 


dal 





pet, PROBLEM FORMULATION 


inviscid adiabatic flow of a gas along the positive x- 
axis is considered which is perturbed by either A) an air- 
meme cascade executing small pitch oscillations (Fig. 1) or 
B) small amplitude vibrations of the surface of a circular 
cylinder whose axis is aligned with the free-stream direction 
img. 2). 

The following basic equations hold: 

The equation of state for a perfect gas, 
> fo) Ral (1) 


the conservation equation for mass, 


Ramee + oy OM) = (2) 
Dt OX dy og 


the conservation of momentum 


poe ee = 16 
Dt for Bios 
Ye ee 
DC p oy 
Oe har 
phe 0 32 


and the conservation of energy 


Ds 
Dt 


= 0 


D 
where Dt is the substantial derivative and s is the entropy 


of the fluid particle. 
The assumption of small amplitude vibrations makes it 
possible to introduce the small perturbation concept. 


aRZ 








Cascade Configuration 





Pigure 2. 


Cylinder Geometry 
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Therefore, all flow guantities are assumed as small pertur- 


bations linearly superimposed on free-stream quantities, 


1.e., 
ee tu (5) 
V=zv (6) 
W = w (7) 
@=c, +c (8) 


and the pressure and density perturbations are: 

Ap = P-Po (9) 

a — 0 cy (10) 
The pressure perturbation can be expressed by the velocity 
Semecound perturbation, l.e., 

P"Po = — PoC (C-Cep) (11) 

Specializing to planar flow, one has for the continuity 
equation: 


eee sc, OM + ce, 2 = 0 (12) 


v-L ot y-1 ax aX oY 


ma the Euler equations: 


eee ee, OS = (13a) 
at aX fate ox 
Beet Sc, OF = 0 (13b) 
at ox yr1 oy 


Equations (12) and (13) are the governing equations for the 
case of supersonic flow past oscillating cascades. 
For the case of supersonic flow past vibrating cylindrical 


shells, the equivalent single equation for the perturbation 


14 





M@eeential, ¢(x,r,6,t) in cylindrical coordinates was taken 
as the more convenient basic governing equation, 1.e., 
cl I 1. 


(Le ase nee ee 


pa cia y 3 OO om zt ce: aac (14) 


eee 


oO 


For a derivation of this equation reference is made to one of 
the standard texts, e.g., Bisplinghof, Ashley, and Halfman, 
"Aeroelasticity", (Ref. 1). 

For a complete formulation, proper initial and boundary 
conditions must be prescribed. The most important ones are 
A) the flow on the body surface must be tangent at any instant 
in time and B) Sommerfeld's radiation condition must be satis- 
Fied, 1i.e€., waves must propagate away from sources of 
disturbance. 

After properly linearizing the general flow tangency 
condition, (see Ref. 1) one obtains for the flow tangency 
condition on a thin airfoil, performing a time dependent small 
amplitude motion about the y=0 plane, 


v= 2hyy 2h 


rae at y=0 (15) 


where h(x,t) is the equation of the airfoil chord line. 
For the vibrating cylindrical shell the flow tangency 


condition over the flexible portion becomes 


el a2 ay, 2B 
oe aX O<x<l 
o| = 0 x<0 (16) 
r=R 


LS 





where h(x,@,t) in this case describes the radial shell 
deflection. 
With the assumption of a perturbation potential, the 


momentum equations, equations (3), may be integrated to 


| 
BeBe tu ont [Poo (t) (17) 
dt 2 ox /?/— 


assuming uniform parallel flow far upstream, 


2 
c(t) ste (18) 
2 


Also, the density-pressure relationship can be expanded as 
i dp an I dp ee tes 2 
p (pb) =p, + ta? (P Po) +3 (G52) (p-p,,) 2+... (1.9) 


therefore, 


p p = rs, 
| se = | p-PBs apeBs oo 
” eo (p) Pp PL py ob 


Hence, the perturbation pressure can be expressed in terms 


OE the perturbation potential, i.e., - 


p=p_= -p, J+ (21) 
dt OX 


In pressure coefficient form this becomes 


Gg cae = 25 - 


: 2 = 4, (22) 


} 
x 


oO 





8 NM 
cG|N 


l 

7. Ue 
For further analysis it is convenient to use non-dimen- 

Sional quantities. Introducing for the oscillating cascade 


case the blade chord as the unit of length one has the 


following non-dimensional coordinates: 


16 





mi 
ll 
Q | ™ 
he 
I 
Q>hK< 


, b= EUe (23) 
eS 
Furthermore, introducing the assumption of harmonic time- 
@ependence, one can define a reduced frequency 
k= — (24) 


where » is the circular frequency of blade oscillation. 
Following previous work by Teipel (Ref. 3) for single 


blade oscillations the non-dimensional complex amplitude 





functions 
Ueoies = * NS) 
Seams} : 
V(z, Fer PF = ——— © (26) 
/M2— co 
Sees, 2 lL. Sc 
C(x,y)e = vet M2 Ce (27) 


are introduced. Substituting these equations into the continu- 


Meyeeand Euler equations (12) and (13) gives 


ac 


eee et OS + ikM2c = 0 (28) 
Ox fe) ox 

oU aC : Es, 

aa te 3x + 1i1kU = 0 (29) 
2 a os ac + ikV = 0 (30) 
x" ygz=T PY 


where the overbar on x and y is being dropped for convenience 
in the subsequent analysis. 

Introducing similar non-dimensional variables with dis- 
tances referred to flexible cylinder length, velocities to 
free-stream speed and time to cylinder length divided by 
free-stream speed, one obtains the equation for.the non-dimen- 


Seonal perturbation potential, i.e., 
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2 i One _ 2 - 
(1-M aa a SS, 4, 6 2M oan se ee = 0 (oa) 
subject to the boundary condition 
o_| = 23 + a eect 
Ee 7 (32) 
= 0 x<0 


If the radial shell deflection is assumed as 


ikt 


h(x,6@,t) = Z(x) cos n6é-e 1 58)) 
where k is now referred to cylinder length, i.e., 
_ woL 
k = ae (34) 
Then the velocity potential must also have the form, 
Siew i,t ,t) = ¢(x,r) cos nefet (35) 


Substitution of equation (35) into equation (31) then gives 


T n2 ; 
(1-M 1 he + 7 + ate za 2ikM b. + k“M<¢ 0* 36) 
subject to 
_ 82 ; 
¢_| = t Tes Dass 1) (37) 
r=R =" 6 
i x<0 


The pressure coefficient in terms of the non-dimensional 
melocity potential is 


oe eee cuako to : So) 
and the pressure coefficient amplitude is 


Cp = -2[ik¢g + d] (38a) 


These are the basic equations used in the subsequent cylind- 


rical shell vibration analysis. 
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ibaa iit tOD OF CHARACTERISTICS 


Pec YLENDRICAL SHELL 

In Section II, equation (36) was developed as the partial 
G@atfterential equation for $(x,r) that requires solution 
Soyect tO the boundary condition, equation (37). The follow- 
ing developments are based on a previous paper by Platzer, 
Brix, and Webster, (Ref. 7). 

If M>1l, equation (30) is hyperbolic and has real character- 


istics, which satisfy the ordinary differential equation 


(M2-1)dr*-dx* = 0 (39) 
An arbitrary function F(x,r) has the following derivatives 


along the two characteristics: 


dF _ dx dr _YM2-1 il 





where ds is the differential arc-length along the character- 
istic ds* = M@dr?2 and 


ax /M7—-1 





ace on er) 
and | : 

she Va 1 

AS = + M (42) 


Letting S,; be the arc-length of the left-running Mach line 
and S$» the arc -length of the right-running Mach line, one 


obtains for the first derivatives, 


YM2=] 1 
M 


Fy) = Ea, (43) 





9 





and 





Fo =, Pa — | (44) 


and for the second derivatives 


ak a 


~-F — pas = (= 
eS Sy TNE [TOTES Boe “Foy (45) 


Therefore, equation (30) can be written as 


d $y 2 Sa iLkM 2_ n2 
mo 0 0 o@eM Ot -j ot Cee KO - 77) ¢ (46) 


| LOOe %, 
NN 
\ 


meeGewonwe Cylindrical Shell Characteristics 
Network 


The finite difference form of this equation for the points 
P and S (see Figure 3) on the right-running Mach line is 


ikM : 


S -_ - i 
Reg = Set gear (eaten) + OR? SFee (47) 


Similarly, for the points Q and S$ on the left-running Mach 


line one obtains 





LkM 
d2(S) - $510) _ $)-¢. _ —— > n2 
AS ma uizan 91 too? + CK owe (48) 
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Now solving for the values of ¢, and ¢5 at S from the known 
values at P and Q (Figure 3) the entire flow field may be 
computed point by point. For this purpose the values for $6, 


and ¢59 in equation (47) are averaged between P and S, i.e., 


o1 = $1 (P) + $1 (8) (49) 


op oad (P) + $9 (S) 
< 2 (50) 


and 
o(S) = $(P) + 5 [o2(P) + $2(S)]As (51) 
is obtained by integration along the right-running Mach line. 
Similarly, $$, and $¢5 in equation (48) are averaged be- 


tween Q and S 


sy PGS me SEL (52) 
Agee $210) + o2 (8) : (53) 

and 
o(S) = 9(0) + 2 [61 (0) + 41 (S)] ds (54) 


This results in a system of two equations for $¢,(S) and ¢$2(S), 


me. , 


Hl 
Q 


Ao 1 (S) + Boo (S) (55) 


Do, (S) + Eds (S) 


ty 


(56) 
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where 


Pes ees bs ikMAs aye) 
4Mr (S) 2 /MZo1 


As evi esac) 
B= -————— ae - ———__—_) (58) 
Mr (S) 2V/M2—-1 M2r2(S) 
As ikMAs 
Caan 1) (1 + —-— 
4Mr (P) o/Me-1 
As ikMAs , As2 n? 


— ba ; = ae ra 
3 do (P) I Atir (P) Dee 1 4 M2y2 a 





(59) 
pe Je eee is) AS? (2 _ Re \} (60) 
4Mr (S) Dal 4 M2r2 (Ss) 
a _ As _ ikMAs 
Be >) 9 Gar (8) ~ Sez Hed) 
As ikMAS 
F = $5(Q)[-1 + + —- 
amr (Q) "9 fZ=T 
As ikMAsS As? n2 
ono — + = - -—(k*? - ——————) ] 
4r(Q)M © 9 vyqzzq 4 M2r2(S) 
Z 
+ 9(Q)[- 52(k? - —B_) - SS(x? - —*—_)] 
r*(S)M2 Mer (0) 
(G2) 


The above procedure will suffice for the computation of 
¢ and its derivatives at a general flow field point, but needs 
to be modified to obtain the values on the initial Mach line 
and along the cylinder surface. 

On the initial Mach line (Figure 3) it is known that 


¢=9,=0. Hence, equation (56) becomes 


E¢o(S) = F (63) 
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where E is known from equation (61) and F simplifies to 


F = ¢2(Q) [-1 + giSay + nas (64) 


On the cylinder surface the normal velocity is prescribed, 


ieee. , 


02 


Pee (Oo) = yah Ike (65) 


be 


ho Ss 


whereas eguation (55) must hold along the right-running Mach 


line. The values for ¢,; and $5 on the surface are then found 


from 
C+ = (2 + ikZ)B 
61(S) = i. (66) 
and 
C - a (32 + ikzZ)B 
62(8) = ot (67) 


The pressure distribution over the vibrating cylindrical 
shell is obtained from the linearized pressure amplitude, 
equation (38a), and the generalized aerodynamic force One is 
found from 


re, ¥ (x) dx (68) 


= 
Cnr 3 


where 


Z(x) = 2 A_ ¥_ (x) (69) 
mm om mM 


and Com (*) is the aerodynamic pressure amplitude due to the 
deflection ¥ a Oe) 
These equations were programmed in Fortran IV. The 


program description appears in Section IV-A. 
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B. CASCADE 

In Section II equations (28) - (30) were developed as the 
Memerning system of partial differential equations for the 
Meiysis of the cascade problem. A solution of these equations 
can be obtained using the method of characteristics. The 
following development of this method is largely based on the 
previous work by Teipel (Ref. 3) and Chalkley (Ref. 9), where 
further details can be found. 

Transforming the system into a new coordinate system, 
E(x,y) and n(x,y), and considering now a system of three simul- 


taneous equations in unknowns U V and C,, the equations 


Se c 
may be written 

2 ne = COlemen v1 

EU acOt a Eye + M EC e U,ty n y 
- M*C ny - ikM2C (70) 
Ee + EC e = = Uae - eae ~ ikU (71) 

V+ = = es _ 3: 
oa E tan oa Eye Rule tan a ou Ti es 2) 


If Ups Vee and C, are made indeterminate across lines of 
n = constant, the determinant of the matrix of coefficients 


of the above equations must equal zero, or 
3 Pee? pO a2 
oe Se? £2) = 0 (73) 
The solutions to this equation give the three character- 


#Stic directions in the physical (x,y) plane: 


dy _ 
aa 0 (74) 
d 

7 ie tan a (75) 
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immoraer tO Lind the compatibility relations of the first 
Gerivatives of U, V, and C along the characteristics, equations 
fue =— (7/2) are written along lines of —€ = constant, jn = 


constant, and along streamlines, that is, for 








str(x,y) = constant, oY = 0 (76) 

wea) = Constant, a2 = : (77) 
ey Moa 

mio, ye = Constant, ot = = 1 (78) 
x YM2z=1 


Tf an arbitrary function of x and y in the equations is de- 


motea by £ then, 3 holding str constant becomes 





= oes VED 
Likewise, eae holding € constant becomes, 
of of i of 
(—) = — + = (80) 
Ooerc OX cea Jy 


and, = holding n constant becomes 


) a a (81) 


Combining these equations and rewriting, 


of 


3x (Sx) str oe 
oer (oe oe 
eee ot (2) ) (83) 
Ope | = of eh 
sy 7 LYAL GS), - (84) 


Substituting these relations into the equations of motion 


gives 
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1, ,3U 3U ee. av, _ ,oV 
S l(a e it ina ad es ules Jb), (a) Sree 
+S M2((35), + (25) ] + ixm?c = 0 (85) 
fe, 3U 1. OC 3C ee ae 
a Bas E + aa ae ate a Mae 1 ‘cee EI + ikU = 0 (86) 


(<=), + (G2) 1 + ¢ UGE), - (GE), + ixv = 0 (87) 
Subtracting equation (86) from equation (85) gives an equation 
that is first added to equation (87) and then subtracted from 
equation (87). The result is two compatibility relations that 
must be satisfied along lines of constant &€ and lines of 


Genstant n, 1.@., 








OV aC : 1 

(ox) a (se + ik[V + eae (M2C-U)] = 0 (88) 
oV oC c« o 1 DZ ae = = 

fee a, tT tklV it ie) a 0 (89) 


The third compatibility relation is obtained by considering 
equation (86) along a streamline 


au 
ox 


ac 
Ox 


(=) es) ae Ca (90) 


str str 
If irrotational flow is assumed, the Sqlatlonwor Irrota— 
tionality 


I (91) 


dy dx 
may be substituted for the second equation of motion in the 
System. With the assumption of simple harmonic motion, it 


May be written in terms of Teipel's amplitude functions, 


oS _ ee ee (92) 
y yx 


The non-dimensional system of equations can then be written 


as, 
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Uva Me = + GkM2c = 0 (93) 
OX Jy ox 
Ue den. _ 
ey MZo oO (95) 
ay ax 

meisanstorming equations (93) - (95) to the arbitrary co- 


ordinate system, & and n, the system can be considered one of 


ES and C.- By making any of 


these derivatives indeterminate, the results of setting the 


simultaneous equations in Ue, V 


determinant of the matrix of their coefficients equal to zero 
are equations (74) and (75), the characteristic directions. 
By writing the equations of motion, as before, along 
lines of € = constant and n = constant and performing the 
same algebraic manipulations, the two compatibility relations 


can be determined: 








aU aV M2 > 
ie 
: 2 
() + (2) + ik —— (u-c) = 0 (97) 
x dx 2 
Mé~l 
The third compatibility relation is, as before, 
9U aC Ree 0 
(oe) stere oe sto (90) 


The equations of motion are now reduced to a system 
containing only derivatives with respect to x. It is, thus, 
possible to solve these equations using finite differences in 
the following manner: separating the real from the imaginary 


parts, the equations can be written, 


Ze) 

















R R 7 
( ores : ( Jae KUL a (98) 
aU L , oC, 
(Sx ) ae (< Mee 33 KU, = 8 (99) 
3U 3V ; 
R R M : 
or? E an E k ; (U_. Cr) = 0 (100) 
Me-l 
3U 3V 5 
i i M _ 
(sx) £ ia 5 1s es (U, Cp) = 0 (Om) 
oU aV 2 
R R M - 
a a ie serra i k on (UL. Cc.) = 0 (102) 
3U 3V 5 
1 I M 7 
a y + Vinca +k oe (U, C,) = «0 (103) 


To write these equations in finite difference form the 
computational molecule shown in Figure 4 will be used. With 
the velocity values all known at Pj), P 5, and P5,, all the 


values at Poo can be calculated. The distance Ax is arbitrary. 





x 


Figure 4. General Flow Field Molecule 
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Tf F is an arbitrary flow quantity, its partial deriv- 


atives can be written as 


ee i 
ax Str Ax (104) 
| Fo, - F 
QE) = 227 Fa on 
1 ax 
va 
Fi, - F 
(22) = 442) | (106) 
2} = AX 


2 


The values of the flow quantities in the equations are 


averages, such that 


fi = > tFy5++ F.,) (108) 
(F) = = (Fo, + Foo) (109) 


With these relations substituted into equations (98) - (103), 


the system becomes after some algebraic manipulation, 


eo opt ~ Kz (110) 
eee Torn ~ X27 (111) 
pee ee 1597 ~~ Yoor) = Kap (112) 
eon ot 950R | C20R) — Saar (113) 
eeeeeo2n) 1 '°22T ~ Voor) = Kser (114) 
mee oe) 7 U5oR ~ So2R’ = Keer (115) 
where 
“oe ~ “abe 4) SU a ne (116) 
“spar Sh sic <0 Saini emia ae (117) 
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K- = 0) - V + BL (U - C ) (118) 





parigeoR 7 “ioR ere oT 

eam ies 9 ‘i2t “1 '°12R ~~ “12R? (119) 
Toran 2am? Yor 7 “217° (120) 

ome 2ir | ote 9 1 21R 7 ©21R? (121) 

and 

A, = > kax (122) 

B. = ik s Ax Ie) 

M2-] 


As shown in Teipel (Ref. 3), solving these equations 


gives, 
} 
me 1-28) (5 (K, tke og) -ByKy 971423, (5 (Ky pth gy) tBLK) op! 
= 2 2 
(1-A,B,)? + (2B) (124) 
_ 7 L L 
Usaz = (1A BL) ID (K3 47 te 67) +B Ky op] - 2B, le (Ka yp Ko gp) By Ky oz) 
es Z Z 
(1-A,B.)? + (2B,) (125) 
Le = 
Vo2R — (KEER K34R) (126) 
aoa 2 
rma Acie “547 (127) 
Coon = Kiar 7 Voor + AUaarx (128) 
eee = AE (129) 


Zu VW Aa Pfeil IT°22R 


These are the finite difference equations for a point in the 
general flow field. 
At a point on the upper surface of a cascade blade, the 


normal velocity Vo is prescribed by the movement of the blade. 
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This velocity was given in the problem formulation. Here the 


Semoutational molecule is shown in Figure 5. 


yA 





Ber J Poo 


Beotweso. Upper surface Molecule 


The applicable equations are, 


ee eRe 201 = “12K oun 
Wage oie Vr) anid era) 
ore oor 1-227  “s6Rr  “34R (te) 
een © 1 °22R «Ser C341 (2S) 


where Kio and Keg are aS before, but, 


Kear — Voor (134) 


Kyat = Voor (135) 


as given by the boundary conditions. Solving for the flow 


quantities gives 


U = (1-A,B,) [K K 


22R ]+2B,[K 


56R §34R PrSiar Be SEAES at Le) 
Za 2 Z 
(1 A,B) + (2B) (136) 
U -« (l- = _ 2 a 
epee Py) KE 6 WK3 47 4B 7% 21-28, [Kae pKa yp 7 BK 07! 
iS 2 2 
(1 A,B.) 4 (2B) Clis7 ) 


C = K =U f: A,U (AES) 


EXINRS JARS ek 220 


Sal: 





& = K ee ey EN roo) 


Ta 
and Voor and Vooz are known from the boundary conditions. 
These are the finite difference equations for the flow 
velocities on the upper surface of a cascade blade. 
Likewise, at the lower surface of a cascade blade, the 
normal velocity, Vo5 is prescribed by the movement of the 
blade, and was given in the problem formulation. The compu- 


tational molecule, different from that at the upper SUR Lace 


of the blade, is shown in Figure 6. 





x 
Figure 6. Lower Surface Molecule 

Here, the applicable equations are, : 
pruimee on  “t°221 ~ “i2R ae 
eT 1 422R ~ “127 ae 
ee aoe 1-291 5 sqR 6 “34R ‘1-7 
ere Doe T-20R = “sqi * “say ‘14> 

where Ki and Key are aS originally stated, but, 
-SGREn 22 aa 
“gat gy \ojohe ait 
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asegiven by the boundary conditions. Thus, the finite 
difference equations for the flow velocities at the lower 


Sumrace Of a cascade blade are, 


| Beene Sort cine ear! e561 * 8347 PX oR] 
a (1-A,B,)? + cae (146) 
ia Stee Gott 12k! <1! SoR 34k oT 127. 
aie (1-A,B,)? + (2B,)? (147) 
ORMIOR ME 22R 10227 (ese) 
2 22T 22k sae 
where Voor and Voorz are given from the boundary condition. 


To obtain the conditions on the initial left-running 


Mach line, assume that P and P as shown in Figure 4 are 


a 1: Zi 


Just in the free-stream and then let Ax shrink to zero. 


Since Ky 5 = Kee = 0 and AL and I go to zero as Ax goes to 


eee, che initial finite difference equations for §=constant 


are, 
U,> + C45 = Q C501) 
U>> - a = Oo . (151) 
U,> + Vo55 = Q (SZ) 
thus, 
ee = 92 «C22 oe 


Equation (96) can be written, 


(), ss ikU (154) 
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Integrating gives, 














M2 
-ik ? *: = 
eal ye ge) = Mem1 (155) 
eauSs , 
eee (0) (ae eee ev ee(O)sin (k x) (156) 
U = = cos 36) Ne Sin Se 
22R 22R ea AG aa 
U V (0) (k M? )+V (0)s3 (k M? (OLS 70) 
=— COs SC ya Sin Da 
DOT DoT im 22R ae 
Voor ~ ~ Yoor (158) 
Vos = = Ue (159) 
Coon ~ ~ Uoor (160) 
ne ee ear (161) 


For the initial right-running Mach line, Pia and Pio Cr 


Figure 4 are assumed to be just in the free-stream and Ax 


ms brought to zero. Since Ki > = Kx, = 0, and A, = ES = Q, 
Uny + Con = 0 (162) 
Us5 - Von = 0 (163) 
Un ~ Un, + Von = 0 (164) 
thus, | 
wee O27 «(22 iro! 


Equation (97) can then be written, 

SU 

Some) an 1LkU (166) 
Integrating gives, 
M2 


M2-1 


IEG 





U = v(0)e ‘i (167) 
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Thus , 


M2 














Usor = Vo op 0) cos (k me x)+V55, (0) sin (k aie x) (168) 

ie = V2.,(0)cos (k = x)-Vo5p(0)sin (k ew (G5) 
M2-1] M2-1 

oe Soe Lee 

oR) ane 22 (oH) 

Von eer ae 

pe MoD (nu) 


At the first position encountered on a blade other than 
the first blade there are two grid points for which the flow 
field quantities must be computed; one on the upper surface 
and one on the lower surface of the blade. 

As in Figure 7 the blade is assumed to protrude a very 


small distance past P thus necessitating the computation 


De 


of the upper and lower flow velocities at this point. In this 


case, the equations developed earlier for an upper surface 


grid point are used for P and those developed for a lower 


22U 


SurLace grid point are used at Poor, 





12 Blade 


Figure 7. New Blade Molecule 


oo, 





Mea pOSition just aft of a blade there are also two 
memaepointsS to be accounted for; one above and one below the 
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assumed to be a very small distance aft of the blade's trail- 


Seo plane in the wake. As in Figure 8 the point P 


ing edge. 


Wake 





Wake Molecule 


Figure 8, 


At this position the equations for an upper blade grid 


Bomnt apply to P and those for a lower grid point apvly to 


22U 
Poor! 1.e., 
or. were ~ se22tu ~~ “12R0 eT 
emer mon | “12TU oN 
Pepe  oiu ee 22tu = “ser “22nu ‘276 
eee one 22kU) ~ “ep, ~ Yooty (477) 
and 
‘logen “° Copten, Sloan, mumps ar Sats) 
oa” “phen | “ocop bre, shee iota 
Usorn * Br Coats, 7 Yoars) = Voorr + Kzapn (189) 
eee ons “Joan, “20rn 7 “34x ‘“*8)) 
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where 


Kiornu ~ Yiiru * Cirru * Ariaiszu (182) 
me ony “lity “1 11RU (183) 
Primemeovirin ts iiRt, © “5 -11TT, (184) 
MUM Gt. | 111 “tT 1IRL (185) 
Veena aor” yi on 7 127? (186) 
K3qr = Uyor ~ Vyor 7 By Uiop ~ Crop) (187) 
Poem “oir * Py \Yoq7 ~ Coir? (188) 
Kogr = Your + Varr ~ Bz oan ~ C2iR? (189) 


Equations (174) (181) form a system of eight equations with 
twelve unknowns. However, since it is known that the wake 
cannot support a pressure difference and that the normal 


velocities on either side of the wake must be equivalent, i.e., 


[Conus (22RD = “20R ee D) 
Come 225L 7 22% | 7) 
Cooru = C22rn = ©22R (oo) 
Smt 0 T ~297 sae 


we can eliminate four of the unknowns and achieve a determinate 


system of equations. Equations (174) - (181) become 
ezorumeeooRes “r4ooTu -— *iorv Gh2g4 
e2enun 4s - 207 * “s95>Ry = “i2tu (195) 
peu Dee 7 “991 ) Yoory’ = Kser eo! 
oe eee t ooRG ~ Coon’ = *se6x oy! 


Shi) 





oreo) 1 227 1 2RL (198) 
Dm 22RL (1271 (199) 
Urort ~ Yoon + By (Soar ~ Voorn) = Kaa (200) 
Perm 220° 1 !22Rt © “22R’ ~ S3ar (201) 


The pressure distributions along the surfaces of the blades 
can be determined from equation (11). Dividing both sides of 
this equation by p_ and using Teipel's amplitude FUNCELONE fem 


c, the equation can be written, 


ies 2 ae = 
Ie «,y) = YM C (x,y) eor ; y=0 (202) 
where 


P(x,y)ero* ro 8 (203) 





Thus the non-dimensional pressure distributions along the 
upper and lower surfaces of the blades in a cascade can be 
Semouced in a point by point calculation of U, V, and C along 
Memes Of €=constant and n=constant. ‘ 

These equations were programmed in Fortran IV. The pro- 


gram description appears in Section IV-B. 


C. THE LINEARIZED METHOD OF CHARACTERISTICS OF SAUER AND HEINZ 
Sauer and Heinz (Ref. 11) developed a linearized method 

of characteristics for supersonic flow past bodies of revo- 

lution at zero angle of attack. This method also is derived 

from the linearized potential equation for steady supersonic 


muOwW, 1.e€., 


CUS) (ya in (204) 


rr 


Hye 
£}- 
i 
© 





Sauer and Heinz showed that the compatibility relations can 


be written as 











Reape yea 7 OU 
SE YM‘-l xr aE (205) 
and 
a(rv) _ yqz=T yr SU (206) 
on on 
where 
_ od 
UP ee (207) 
_ ao 
sip a (208) 
and 
FE=x - YM¢-lr (209) 
eee VMe—) 7 (210) 


are the characteristic coordinates as explained previously. 
Petting this equation in finite difference form one 
obtains for u and v at point 4 (Figure 9) using the known 


values in points 2 and 3, - 





(rv) ,-(rv) .t+V7M2-l(r utr u 
3 ee (211) 








4 
2.2 
YM2—-1 ea 
r.(rv).tr (rv) ,.-VYM@2-1 rir (u.-u.) 
(vr), = 2 sale 2 E23 (212) 
r +r 
E n 


where oe and r, are the average r coordinates between 3 and 
4 and between 2 and 4 respectively. Computing the flow pro- 
perties at the body using the exact flow tangency condition 


leads to 
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dh 





YM*-1] r u,t+(rv),-r- (=) 
Guin 
Pe uaa 
yYM@-l or + Us (ay) 5 
_ sb | 


where r = h(x) describes the body contour. 


However, in order to apply the Sauer-Heinz procedure to 


the case of an infinitely long cylinder with a wavy surface 


between 0O<x<l a different starting procedure is needed. 


is eaSily seen from Figure 9 that 


h aL 
ey) a0) ( ———_—_—__—_—_—_—_—_—- + ]) (215) 
- 2 YM7-1 tan So 
realizing that 
R. + R 
eet! 2 
i = (216) 
where 
R, = Ro Fie Ze) 
and 
R. = Ro + hy (218) 
one obtains 
Ao i 
r =R+t We (3+ ———___.._____._ } (219) 
YM7-1 tan 6 


0 
Since there is no disturbance upstream“of the leading edge. 


Mach line, one obtains from equation (213) 


~ (Rho) tan 5 9 
0, = —_——_————————— (220) 
YM2-1 rit (R+ho) tan § 5 


Pmere the point 2 is the first body point. 
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Finally, the pressure coefficient is obtained from 


= 2u - v2 + uz (M2-1) fin) 


The program listing for these equations is given in 


Appendix I. 





As, 
ast 
R 
ss Oe 
x 
Figure 9. Sauer-Heinz Characteristic Network 
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Lae COMPUTATIONAL PROCEDURE 


fey LL NDRICAL SHELL COMPUTER PROGRAM 

This program calculates the pressure distribution over a 
memmonically vibrating cylindrical shell. It is valid only 
for a cylinder in a supersonic flow field vibrating with 
small amplitude vibrations. 

The program calculates the flow field properties adjacent 
memine Cylinder using the method of characteristics finite 
difference equations developed previously in Section III. 

The pressure distribution is then integrated using the trape- 
ZOidal rule according to equation (68) to obtain the general- 
ized aerodynamic force. 

The computational diagram shown in Figure 3 is used to 
calculate the flow field quantities at S in the general flow 
field and on the cylinder surface. The distance As shown in 


the figure is determined by 


pescosse (222) 


where v is the grid fineness ratio. The grid fineness ratio 
is an arbitrary input variable into the program and is equal 
to one less than the total number of grid points on the 
cylinder surface. 

Input parameters are entered into the program on two data 
Seas. The first card contains the date, e.g., 11 APR 1973, 
in the first twelve columns. The second card contains the 


following input data: 
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oeaCk DATA 





1-20 Free-stream Mach number, M 

Pa —30 Radius to length ratio, R 
31-40 Reduced frequency, k 

41-45 Axial mode number, m 

46-50 r, as defined in equation (68) 


for the generalized aero- 
dynamic force. 


51-55 Circumferential mode number, n 
56-60 Grid fineness ratio, v 
where the first three quantities are floating point numbers 
and the rest are integers. The definitions of these para- 
meters are given again in the alphabetical listing of all 
program variables in Appendix A. 
The program has five subroutines. They are as follows: 


COMPXR Computes the x and r position 
of the given grid point. 

MACHLN Computes the flow quantities 
along the initial left-running 
Mach line at the given grid 
joleniiguer. 


GENFPT ' Computes the flow quantities at 
a general flow field grid point. 


PANPNT | Sctpuces Ene Flow quantities on 
the cylinder surface at the given 
Gia DOa nie. 


INTEG Computes the generalized aero- 
dynamic force. 


Flow quantities along the initial left-running Mach line 
are computed using the equations derived in the method of 
GMaracteristics for initial conditions. Those quantities on 


the cylinder surface are computed using the equations derived 
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meet ne boundary condition. The values of the velocity 
Meemiital and its derivatives are derived as a result of the 
analysis of the flow tangency condition and have been given 
previously. The procedure starts at the initial Mach line 

and continues down succesSive right-running Mach lines. The 
pressure coefficient is computed when the cylinder surface 

is met and stored in an array for later use in the integration 
of the generalized aerodynamic force. 


A flow diagram of this program is given in Appendix B. 


B. OSCILLATING CASCADE PROGRAM 

This program calculates the pressure distribution over 
the upper and lower surfaces of the blades of a subsonic 
leading edge locus cascade which is oscillating in pitch at 
arbitrary frequency and interblade phase angle. It is valid 
only for a cascade in a supersonic flow field vibrating with 
small amplitude vibrations. The program calculates the flow 
field quantities within the area governed by the aft-running 
Mach lines emanating from the leading edge of the first blade 
and the forward-running Mach lines emanating from the trailing 
ee] Of the last blade of the eee (see Figure 10). These 
flow field properties are calculated using the method of 
characteristics finite difference equations developed previously 
in Section III-B. 

The computational molecules shown in Figures 4, 5, 6, 7, 
and 8 are used to calculate the flow field properties in the 
general flow field, upper surfaces, lower surfaces, leading 


Begs Of a Staggered blade, and the wake aft of a blade, 
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/eepectively. The distance 4x shown in the figures is 
determined by 


AX = 2S cot a C22 


fee is the grid fineness ratio. This parameter is an 
arbitrary input variable into the program and is equal to one 
less than the total number of grid points on any Mach line 


bet vans "from one blade to the next. 


ss MD 
\ 


“TS 


Pectmem!O., Cascade Computational Area 


Input parameters are entered into the program by means of 
the NAMELIST option of Fortran. In using this option, input 
parameters are combined in a titled list and then referred 
to in the program by this title (NAM1). The real advantage 
Of USing NAMELIST is that each data card need only contain 
the value of one input variable while the format of the data 
card is such that the variable's name is punched on the card 
as well as its value. In this way, when a different value 


for an input variable is desired, rather than changing the 
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entire data deck, only one easily identifiable Sislies| (eau selerate! 
be changed. 

The format for a NAMELIST data deck is as follows: The 
mest column in each card of the data deck is left blank. 

The first card starts in the second column with the symbol & 
(ampersand) followed immediately by the title of the NAMELIST 
(here, NAM1). The next cards list the variables and their 
iee@ue Values, one per card, in any order. The format here, 
starting in the second column and followed by a comma, is: 
variable name = value. The last card in the data deck and the 
Signal to the computer that the NAMELIST has ended again 
starts in the second column with the symbol & immediately 
followed by END. 

The value of an input variable may be written in any 
format as long as integer variables are written without deci- 
mal points and real variables are written with a decimal point. 
Blanks are taken as zeroes, but a comma must appear somewhere 
between a desired input value and the next data card. The 
date of the run is entered on a data card in the first twelve 
Spaces. This data Berd 1s inserted prior to the NAMELIST 
dataset. 

To illustrate the input procedure, the following example 


is proposed: 


we JONG d= 0.8 

M= ¥2 g = 55° 

ro O79 5 Si= 452 
Xo = 0.5 

5 blades 
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The data cards needed are as follows: 





CARD FORMAT , DESeCRLPTLON 
i. POON LOY 3 Date 
2 &NAMIL NAMELIST title 
S NGRDFN=100, Grid fineness ratio, 


must be an even integer, 
Subject eto: the Constrammt 
of equation (227). 


4 2 FSTRMN=1.414214, Free-stream Mach number 

5 REDFRO=0.80, Reduced frequency 

6 XSUBO=0.50, Elastic axis position 

q - TNWDST=0.80, Vertical distance between 
blades. 

8 STGANG=55.0, Complimentary stagger angle, 
bans 28COt a. 

S FAZE=45.0, Interblade phase angle 

nO NUMBLD=5, Number of blades 

ll SEND End of NAMELIST 


The definitions of these parameters are given again in the 
alphabetical listing of all program variables in Appendix C. 
The program has nine subroutines. They are as follows: 


INPUT Reads in all input data and sets up the finite 
difference grid. 


iNT TAT, Initializes all the flow quantities at the leading 
edge of the first blade and initializes most of 
the logic variables. 


MACHLN Computes the flow quantities along the initial left 
and right-running Mach lines at the given grid 
WO le 

TOP Computes the flow quantities on the upper surface 


of a blade at the given grid point. 


47 





GENFPT Gemoutes tne flow quantities at the given general 
Blows Gaea point. 


BOTTOM Computes the flow quantities on the lower surface 
of a blade at the given grid point. 


NEWBLD Computes the flow quantities at the leading edge 
of all staggered blades. 


WAKE Computes the flow quantities in the wake aft of 
each blade at the given grid point. This sub- 
BouLtine wses LEI subroutine SIMO for the solution 
of its simultaneous equations. 


COMPXY Computes the x and y position of the given grid 
point. 


In INPUT the grid is determined using the grid fineness 
ratio, blade spacing, number of blades, and the free-stream 
Mach number. The value of Ax as shown in equation (223) is 
obtained. The stagger angle is then adjusted such that the 
distance the upper blades are staggered back is the nearest 
multiple number of Ax increments. The "compatible stagger 
angle" is then printed out. The blades may be staggered such 
that the first staggered blade is at least one Ax increment 
aft of the initial left-running Mach line, but no further 
merward, Or, 

d é (Eanes = “COL aja. (2a) 

Flow quantities along the initial left and right-running 
Mach lines are computed uSing the equations derived in the 
mem~mod Of characteristics, Section B, for initial conditions. 
Those flow quantities along the blades' surfaces are computed 
using the equations derived for the appropriate boundary 
condition. The equations for the normal velocities are de- 
rived as a result of the analysis of the flow tangency con- 


dition and have been given previously. The flow quantities 
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in the wake aft of each blade are computed using the set of 
simultaneous equations derived as a result of the assumptions 
made about the flow quantities on either side of the wake and 
are given in the “Method of Characteristics, Section III-B. 
The procedure starts at the leading edge of the first 
blade and computes the flow quantities at grid points along 
@@emrnitial right-running Mach line. It then jumps back to 
the initial left-running Mach line and computes the flow 
quantities along the next right-running Mach line, using the 
applicable subroutines at the appropriate grid points. Only 
flow quantities on the Mach line presently being computed 
and the Mach line which preceeded it are kept in storage at 
aoyeone time. Flow quantities at successive grid points are 
computed in like manner until the trailing edge of the last 


blade of the cascade is met. 


A flow diagram of this program is given in Appendix D. 


feo AUER-HEINZ COMPUTER PROGRAM 

This program computes the pressure distribution over a 
Syeindrical shell or Dodvwaer Devolltionem it 1s written so 
that the shape of the body or cylindrical surface may be 
easily enced: It is valid only for non-oscillating bodies 
Or shells in a Supersonic flow field at zero angle of attack. 
The following discussion applies to the cylindrical shell. 

The program calculates the flow field adjacent the 
cylinder using the method of characteristics finite difference 


equations developed previously in Section III-cC. 
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The computational diagram shown in Figure 9 is used to 
@empuce the flow field properties at 4 in the general flow 
field and at 5 on the cylinder surface. The distance As 
shown in the figure is determined by, 


As = <a (225) 


where v is the grid fineness ratio. The grid fineness ratio 
1s an arbitrary input variable into the program and is equal 
to one less than the number of grid points on the initial 
Mach line. The length Ast n=1,2,3... 1s determined by the 
intersection of the right-running Mach lines with the surface. 
Input parameters are entered into the program on one data 


eee This card contains the following input data: 


SPACES DATA 
1-10 Free-stream Mach number, M 
i= 2 0 Maximum amplitude of the axial mode, Ao 
ZA 30 Gvyilaneer radius to length ratio, R. 
31-40 Axial mode number, m 
41-50 Grid fineness ratio, v 


@meoes the first three quantities are floating point numbers 
and the last two are integers. The definitions of these 
quantities are given again in the alphabetical listing of 
all program variables in Appendix E. 
The program has three subroutines. They are as follows: 

COMPXR Computes the x and xr position of the 

grid point and determines the position 

of the intersection of the current 


right-running Mach line with the body. 


GENFPT Computes the flow quantities at a general 
Plover tela grid point. 
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SURFAC Computes the flow quantities at the 
given surface point. 


Flow quantities on the initial left-running Mach line 
are assumed to be zero. Those quantities on the surface are 
computed using the appropriate equations developed in 
Section III-C, for surface grid points. 

The procedure starts at the initial Mach line and con- 
feneres VOW suctessive ‘right=running Mach lines. The pressure 
coefficient is computed when the surface is met. 


A flow diagram of this program is given in Appendix F. 
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Ve RESULTS AND COMPARISONS 


fee CYLINDRICAL SHELL 
Pressure distributions and generalized aerodynamic forces 
were computed for sinusoidal deflection modes, 1.e., 


Ya = sin (mtx) Westie 2, ae. (226) 


Figures 11 through 16 show results for the generalized aero- 
dynamic forces O14 and Q15 as a function of circumferential 
mode number for two radius to length ratios, r=0.5 and r=1.0, 
@tea Mach number of v2 and reduced frequencies of zero and 
one. Also shown'is the comparison with two previous solutions 
by Dowell and Widnall (Ref. 2) and by Matsuzaki and Kobayashi 
(Ref. 5). The Dowell and Widnall solution is also based on 
the complete linearized unsteady potential equation, but uses 
Laplace transform techniques for its solution. As can be seen, 
agreement between this solution and the method of character- 
memes SOlUtion is excellent. For the sake of interest, a 
recent approximate solution by Matsuzaki and Kobayashi, valid 
only for large values of the circumferential mode number is 
also shown. 

Figures 17 through 24 show the generalized aerodynamic 
forces Q 


and QO as a function of radius to length ratio 


Wee IL AAR 
again for a Mach number of 72, zero reduced frequency, and 
circumferential mode numbers n=0,1,2, and 3. The main purpose 
of these figures is to show that the infinite cylinder theory 
1S a poor approximation for small radius to length ratios, but 


approaches the exact theory for large values of R. 
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Finally, Figure 25 shows the change in pressure distri- 
bution as a function of radius to length ratio for the case 
f-v¥2, m=1, n=0, k=0. These distributions were computed using 
the characteristics method outlined in Section III-A and the 
Sauer-Heinz theory as described in Section III-C. Complete 


agreement was found between the two methods. 


B. CASCADE 
Figures 26 through 29 show the instantaneous pressure 

distributions on the first and second blades of a finite 
cascade for four different values of reduced freguencies. 
The pressure difference between lower and upper blade surface 
is plotted for the instant when the airfoil is at its mean 
horizontal position and pitching nose up. Results for these 
cases were given very recently by L. E. Snyder, of Pratt and 
Whitney Aircraft (Ref. 8), for the following parameter 
combinations: 

M = e175 

See = cs etd = 0.317) 

Stagger angle =@5° (gp = 65°) 

Xo = On.5 


S a0? (Blade 2 leads blade 1) 


Lem 04 ,0.33270.632,1.25 
It can be seen that appreciable differences exist ina 
number of cases, in particular in the wave reflection regions. 


These differences are not fully understood at this time. 
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Figure 27. Cascade Pressure Coefficient Difference 


VG uP Gacimken Chermd 


70 





1.0 4 





0,5 
A BLADE 1 Bee 
f 7 
y , ya 
0.0 cs aX 
yp 
=a 5 & a 
i> — 
dy a _ oe 
<a 
SS SSN VlER (ees, 2) 
154 O GARRICK-RUBINOW (REF, 12) 
PRESENT STUDY 
7 BLADE 2 M=1,5 
Sn za RE), 252 
A\c §=l10" 
“7a | 
0,5 | 
| ZA 
Z 
| 
0,0 / Y 
Aes 


0,0 OFZ 0.4 0,6 0.8 nd 
FRACTION CHORD, x 


Figure 28. Cascade Pressure Coefficient Difference 
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Snyder uses a finite difference procedure and further work is 
required to evaluate the validity and accuracy of this pro- 
cedure is compared to the method of characteristics procedure 
developed in the present work. However, a check could be made 
for the first blade solutions for which the method of singular- 
ities of Garrick and Rubinow (Ref. 12) is available. This 
solution was programmed according to the eguations given in 
Appendix H and the results are shown in Figures 26 through 
29. The method of characteristics solution is in complete 
agreement with the Garrick-Rubinow prediction whereas consider- 
able deviations are seen to exist for the finite difference 
ereaictions. 

A further check using the method of Singularities was 
developed for two blades oscillating with zero phase difference. 
This approach is also described in Appendix H and a comparison 


with the method of characteristics is given in Table 1] below: 


TABLE 1 


CE a 
» Ea ee ee ee 
aoe NRoOmecnliNolc. Miolcuminmonce mM.o.s. |M.0.c. M.O.S. 







Oma Oo Oo or Oo” oro = 
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Further method of characteristics results for a finite 
subsonic leading edge cascade are shown in Figures 30 through 


32 for the following parameter combinations: 


M=/72 Xp)=0.5 
Glin) k—@r > 
g=59.8° Sal BIN 


As mentioned previously, the computer program is written 
sufficiently general so that the pressure distribution can 
be computed and printed out for arbitrary specified Mach 
number, reduced frequency, stagger angle, blade spacing, and 
number of blades subject to the following constraint: 


Te 1) (cance cot nes: =1<800 (227) 


Equation (227) is the expression which governs the maximum 
grid fineness that may be used with a given set of input 
parameters. This constraint is necessary because the program 
is dimensioned such that the maximum number of grid points 


allowed on any one right-running Mach line is 400. 
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Nel. CONCLUSIONS AND RECOMMENDATIONS 


A solution of the complete linearized unsteady potential 
equation has been developed for supersonic flow past vibrating 
cylindrical shells using the method of characteristics. 
Pressure distributions and generalized aerodynamic forces 
were computed for shells of arbitrary radius-to-length ratios, 
axial and circumferential wave numbers, Mach number, and re- 
duced frequencies. A solution for this same problem was pre- 
viously obtained by Dowell and Widnall using Laplace transform 
techniques. Excellent agreement was obtained in all cases. 

It is believed that the present method can be applied more 
Saomiy than methods requiring complex variable techniques. 
It should be noted that the procedure can easily be modified 
to include arbitrary axial deflection modes. 

A further solution of the complete linearized unsteady 
potential equation using the method of characteristics was 
developed for supersonic Digvepast yilecaerng tlat plate 
cascades with subsonic leading edge locus. Pressure distri- 
butions were computed for finite cascades of arbitrary 
Stagger and interblade phase angles performing small ampli- 
tude pitch oscillations. Preliminary comparisons with a 
recent finite difference solution by Verdon and Snyder showed 
reasonable agreement. However, further work 1S necessary to 
assess the validity and accuracy of the present method. This 


Should be done by further comparisons with the Verdon and 


he 





Emyaer procedure and Fourier and Laplace transform methods 
currently under development at Stevens Institute of Technology, 
General Motors Corporation (Detroit Diesel, Allison Division), 
and General Electric. Also, flutter calculations should be 
carried out using the method of characteristics predictions 

of the oscillatory aerodynamic forces which would allow com- 


parisons with the Pratt and Whitney cascade flutter tests. 
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AoE ENE A A 


Sooty ORT CAL oHably COMPUTER PROGRAM VARIABLES 


A. LOGIC VARIABLES 
The following variables appear in the logic statements 

Srethe program: 

ITHAVEP The number of the computation along any 
right-running Mach line. IHAVEP is used as 
Ne tibet “SUDSCEIDeE LOmeas Grid pOlnt. 

ISWICH Switch variable used to designate the 
right-running Mach line being calculated, 
alternately set equal to "1" or "2" and 
used as a second subscript. 

JSWTCH Switch variable used to designate the 
previously computed right-running Mach 
line, always the opposite of ISWTCH and 
used aS a second subScript. 

KOUNT The value of IHAVEP at the cylinder surface. 

ML The number of grid points along the initial 
Mach line and along the cylinder surface. 

B. QUANTITY VARIABLES 

The following variables in the program take on the values 


defined below. All are floating point quantities unless 


otherwise stated. The dimensioned variables are listed first. 


AMAGCP (400) The imaginary part of the pressure coefficient. 
CP (400) The complex pressure coefficient. 
mat (400, 2) The derivative of the velocity potential 


along the right-running Mach line. 


PHK (400,2) The derivative of the velocity potential 
along the left-running Mach line. 


QOPHI (400,2) The velocity potential. 
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PBeALCP (400,2) 
pope (400,2) 
oe ;c,D,E,F 


BETA 


DELTAS 
PST STR 


FINGRD 


HDSTRL 
HEIGHT 
ian 

IKM 
IKW 


IPHI1, IPHI2 


K 
MACH 


MANGL 


NF REQ 
NR 


PARW 


Pi 


ex 


RADIUS 


The real part of the pressure coefficient. 
imemoeednad © COerdinates of a grid point. 
See equations (57) - (62) 


The square root of the guantity Mach number 
squared minus one. 


The step size along the characteristics. 
The horizontal distance between grid points. 


Phe grid™iineness ratio. 
less than 400. 


An even integer 


One-half DSTSTR. 
The increment HDSTRL > tan a. 


The square root of minus one. 
M 


M—1 
Tiesauantity,, 1k Sin mix 





Tine LaACcOL, Lk 


The velocity potential as described in 
equations (51) and (54). 


The reduced frequency. 

The free-stream Mach number 
The Mach angle, a. 

The axial mode number. An integer. 
The circumferential mode number. An integer. 
xr, as Gefined in equation (68). An integer. 


The partial derivative of the deflection 
taken with respect to x. 


Se a595 


The partial derivative of the velocity 
potential taken with respect to x. 


The radius to length ratio. 


imewey lander deflection. 
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APPENDIX 3B 


CYLINDRICAL SHELL FLOW DIAGRAM 


Read/write DATE 
Read Input Data 


Define XLENGT, HEIGHT, DELTAS, 
DSTSTR, HDSTRL, ML, 
TRNGLH, MANGL, II, BETA 














ISWTCH= 2 
JSWTCH= Ll 
IHAVEP= 1 
KOUNT= 1 





Seuouce OPHE(1,1), PHK(1,1), X(1,1), 
ie lei, PARW; PHI (ie), 
CP(1,1), REALCP(1), AMAGCP(1) 

Write Input Data 


KOUNT= 2 





104 
CALL COMPXR 
Compuce PARW, W, AKM, 




















IKM, XKW, IKW 







RmAAVEP SEO, | 


THAVEP .EQ.KOUNT 


CALL GENFEPT 


Compute PX, PR 


1000 
LHAVEP .EQ.KOUNT 


CALL MACHLN 





CALL PANPNT 








ITHAVEP= I[THAVEP+1 
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i Zo . 
iSweCd BO. 1 


TSWTCH=1 
JSWTCH=2 
THAVEP=1 









ISWTCH=2 
JSWTCH=1 
THAVEP=1 


KOUNT=KOUNTH1 


(KOUNT-1) .EO.ML 
EF 
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OSCILLATING CASCADE COMPUTER PROGRAM VARIABLES 


A. LOGICAL VARIABLES 
The following variables appear in the logic statements of 


the program: 


CO Set equal to "1" when a blade or wake com- 
putation is made. Set equal to "0" after 
the next general flow field computation 
occurs. 

ITHAVEP The number of the computation along any 


mECit-ruinning Maen tine Starting WLEh Zero 
Ae cea initiml lere-runming Mach Jaine.- The 
Valves THAVEP+Iges the £irst subscript for 

Bae VOVAECIU AT tae y,.. 


ISWTCH Switch variable used to designate the 
right-running Mach line being calculated, 
alternately set equal to "1" or "2" and 
used as a second subscript. 


JLINE The Mach line counter incremented each time 
Gureic-rumiing Machalines rs —COMp.Leted. 
Reset equal to the value of IHAVEP at the 
leading edge of a new blade. 


JSWTCH Switch variable used to designate the 
pEevyLously comeuted raght—-running Mach 
line, always the opposite of ISWTCH and 
used as a second subscript. 


LCOUNT The value of IHAVEP at a surface or wake 
Gist, poOine-, 

LIMIT The maximum value of IHAVEP. 

LIMW The value of LCOUNT at the trailing edge 
of a blade. 

MAXI The value of JLINE on a right-running Mach 
line which intersects a blade's leading 
edge. 

NEWDST The number of Ax increments in a blade chord. 
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NSTPTS The number of AX increments between the 
initial Mach line and the leading edge 
Gumche @earst blade. 


NUM Set equal to one less than the number of 
the blade encountered on a Mach line. 
NUMBLD The number of blades in the cascade. 
NUMOLD Preserves the maximum value of NUM. 
NUMP I NUMOLD plus one. 
OLDL Preserves the value of LCOUNT at the highest 


blade or wake yet encountered. 


Pee QUANTITY VARIABLES 
The following variables in the program take on the values 

defined below. All are floating point quantities unless 

otherwise stated. The dimensioned variables are listed first. 

A(8,8) Matrix of coefficients in wake computation. 

B(8) Right-hand side of equations (194) - (201) 

fee Z2. ,V22R,V221,C22R,C221(400,2). The velocities, real 
and imaginary of general flow field, upper 
blade, and upper wake points. Equivalent 
temtne respective Te1pel amplitude functions, 


Usop elmo ch) Coot: 


Meee U331,V33R,V331,C33R,C331(400,2). Same as above, but 
uSed only at lower blade and lower wake 





Pelnesa 
my (400,2) The x and y coordinates based on chord length. 
AI A, as defined in equation (122). 
Z 
ARG The argument k Pe 
M*-1] 
BI By as defined in equation (123). 
DELTA The interblade phase angle in radians. 
DELTAS Incremental distance between grid points 


along Mach line. 


DSTSTR The increment Ax. 


85 








FAZE 

FMANGL 
FNGDPT 
FSTRMN 


eer... ,G9 


BDSTRL 


ieee... ,KOOT 


KS 


NGRDFN 


REDF RQ 


STAG 


STGR 


STGANG 
S 

TNWDST 
TRNGLH 


U 


V 


VIPANL 


VRPANL 


W 


The interblade phase angle in degrees. 
Mere ilrcl lice 7). iC 

Floating point NGRDFN. 

The free-stream Mach number. 
Intermediate factors used in velocity 
computations in GENFPT, TOP, BOTTOM, 
and NEWBLD. 


The ancrement Ax/2. 


Kyopre ss Keer as defined in equations (116) - 
C7 ie 


Dummy variable used in SIMQ. An integer. 


The grid fineness ratio. An even integer, 
less than 400. 


The reduced frequency. 


The distance the second blade is staggered 
back from the y-axis: 


The distance from the initial Mach line to 
the second blade's leading edge. 


The stagger angle in degrees. 
The factor YM¢-1. 
The vertical distance between adjacent blades. 


. ' 
AWeve increment ZAXx COS a. 
M2 
M2~1 
M2 
M@-] 





Mies LacteOr “COs (kk Jie 


TMeetactcor Sink Dah ae 





The imaginary part of the normal flow 
velocity at the first blade's leading edge. 


The real part of the normal flow velocity 
at the first blade's leading edge. 
M2 
YM2=1] 


The factor 
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XSUBO 


XLNGTH 


XN 


ANEW 


The elastic axis position based on chord 
Leme eh. 


The length of the initial Mach line between 
the planes of adjacent blades. 


Floating point NUM. 


The x-coordinate on a blade based on the 
blade's leading edge. 
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APPENDIX D 


CASCADE FLOW DIAGRAM 





Read/Write DATE 
CALL INPUT 
OLDL=0 

I CO=0 

CALL INTIAL 
Sseomoute NEWDST, LIMIT, 
LIMW, MAXI 











NUM=0 
NUMOLD=0 
NUMPI=1 


1000 
ALL. COMP CALL MACHLN 
7 
LCOUNT .EQ.0 = LHAVEP .EQ. LIMIT )= 
F 


i 


mA Ee. EO. LCOUNT Z 2 ITHAVEP=ITHAVEP+] 
1S 
Ak 
PAV EP. EO .0 
ley 


THAVEP .EQO.NUMPI*NSTPTS 
- AND. 
JLINE.EQ.MAXI 


CALL NEWBLD 
TCO=1 
NUMOLD=NUM 
LCOUNT=IHAVEP+NGRDFN 
JLINE=IHAVEP 
THAVEP=IHAVEP+1 
OLDL=IHAVEP 

NUMPI=NUM+1 
MAXI=NGRDFN+NUMPI*NSTPTS 
NUM=NUM-1 



















Sebi GENEPT 















PAVE oO. LIMIT 


iiAVEPY. EO. LHAVEP+1 


4 





THAVEP .GE.LIMW 


hy 
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enn bOLrOM 


CALL TOP IHAVEP.EO.LIMIT i" TCO=0 





NUM.GE.NUMBLD-L 
~AND. 
DAVIE EO aig TMi 


F 
THAVEP.EQ.LIMIT }+ 5 
F 


LHAVEP=LTHAVEP+1 


THAVEP.EQ. LIMIT -8 


Jig 
HAY he. EO. LHAVEP+ 


NUM.EO.0 }= 


NUM=NUM-1 
LCOUNT=LCOUNT+NGRDFN 
LIMW=LIMW-NSTPIS 













MOM EO a0 us 
LCOUNT=LCOUNT+NGRDFN 
NUM=NUM-1 





> 
THAVEP=0 





NUM=NUMOLD 
JLINE=JLINE+ 


MowlCH.EO.L . 





ISWTCH=2 
JSWICH=] 
alt 
NUMOLD.GT.0O LCOUNT=OLDL 


Pr LIMW=NUMOLD*NSTPTS+NEWDST 


OLDL=OLDL+1 
COUN. EO. LCOUNT+1 


Sra. 
END 


89 





APPENDIX E 


SAUER-HEINZ COMPUTER PROGRAM VARIABLES 


fm. LOGIC VARIABLES 


The following variables appear in the logic statements 


at athewprogram: 


PAY EP 


ISWTCH 


JSWTCH 


KOUNT 


tie number @iethe computation along any 
right-running Mach line starting with zero 
at the initial left-running Mach line. The 
Verte. SHAVEPTIio 1s the first subscript force 
Pliew “quantity. 


Switch variable used to designate the 
Gught=sumning Mach Line being calculated, 
alrennacely Ser equal to "liom. {2" and 
used as a second subscript. 


Switch variable used to designate the 
previously computed right-running Mach 
line, always the opposite of ISWTCH and 
used as a second subscript. 


The value of IHAVEP at the surface. 


fae QUANTITY VARIABLES 


The following variables in the program take on the values 


defined below. All are floating point quantities unless 


otherwise stated. 
CP (200) 


DELX , DELR (200) 


e200 ,2) 
0200 ,2) 
mk (200, 2) 


A,B 


The dimensioned variables are listed first. 
The pressure coefficient. 


The incremental distances between grid 
points on a right-running characteristic. 


The streamwise velocity perturbation. 
Rhey product, vie. 
PieexanGs COOLTdinaees on thesgrid point, 


Used as factors in the equations for u and v. 
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BETA 


DELTAS 


DELTAX 


DELTAY 


FINGRD 
HO 
MACH 
MANGL 
MFREQ 


NGFN 


NUMPTS 


|eaa 
RETA 
RO 


RZETA 


SLP 
XM 
AMAX 


YO 


The square root of the quantity M¢-1. 


The incremental distance between adjacent 
dmlidpoimes along the anitial Mach line. 


The change in x from one point to the 
next on the initial Mach line. 


Ghee schangesim © from one point to the next 
Pieene wa pViehaleMach lane. 


Floating point NGFN. 

The maximum-v-amplitude of the deflection. 
Free-stream Mach number. 

The Mach angle, a. 

The axial mode number, an integer. 


The grid fineness ratio, an even integer, 
less than 200. 


The number of grid points on the initial 
Mach line. 


Bla Loge 
v6 
n 
ine radsus to length ratio. 
“ie 
tiem renqthn or the aniciral Mach lines 
Slope of the surface at a surface grid point. 
Floating point MFREQ. 


The value of X at a surface grid point. 


The intercept of a right-running Mach line 
With thesge=axis- 


eal 





APPENDIX F 


OSAUER-HEINZ FLOW DIAGRAM 






Read/Write Input Data 
Compute BETA, MANGL, FINGRD, 

NOME Soo bien, 
DELTAX, DELTAY, XM 









ISWTCH=2 
JSWTCH=1 
THAVEP=0 
IT=ITHAVEP+1 






menmouce X(L,1), R(L,1), U(l,1), V(1l,1) 
KOUNT=1 


- D 
foul COMPAR 






THAVEP.EQ.KOUNT 
F 


CALL CENERT 
ITHAVEP=ITHAVEP+ 





U(1,ISWTCH)=0 
V(1l,ISWTCH) =0 
THAVEP=ITHAVEP+1 









THAVEP=0 


fowheCh soe. 









ISWTCH= 
JSWICH=2 


ILSWTCH=2 
SWTECH= 


KOUNT=KOUNT+1 


ie 


CALL SURFAG 
Cao 
F 


5 LOR 
END 








APPENDIX G 


SUPERSONIC FLOW PAST AN INFINITELY LONG 
WAVY WALL CYLINDER 


Leonard and Hedgepeth (Ref. 4) gave a solution for super- 
sonic flow past an infinitely long wavy wall cylinder. De- 
tailed computations uSing this solution were carried out by 
Anderson (Ref. 14). For the cylinder wall deflection 

h = h)sin mmx cos n 0 (Cae) 
where m is the number of half-waves per unit axial length 
and n indicates the number of circumferential waves, one 
obtains for the pressure distribution. 


MW 





Cp = 2ho M Ay Sin(mnx + ¥y) (G-2) 
where 
2 as oo: 
aed _ M He [mnR (M“-1) 2] 
YM7=] 


] 
He 2° fmaR(M2-1) 2] (ey) 


For the actual computation equation (G-3) was rewritten as 





a,et*1 7 ey) oa tee) om 
Dy Chinn On aI | 
(G-4) 
where 
y = (G-5) 
For the case n = 0 the equation was rewritten as 
Ae Joly) - 4 YX, fy) M 
A,e a = = sp ae es (G-6) 
~J, ly) + i Y, (y) YM2=1 


2 ie) 





and 





The program listing appears in Appendix I. 
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(G-7) 


(G-8) 





APPEND AX Ji 


TWO BLADE SOLUTION USING THE 
METHOD OF SINGULARITIES 


As shown by Garrick and Rubinow (Ref. 12), the velocity 


potential of a single blade oscillating in supersonic flow is 





Me 
=1_ =.x-yYM7-1 v(eye Myzay WE) se Ve) SNPS spaiicts 


i, y)= f 
¢ YM“-1 0 M?-1 


(H-1) 





Blade 1 — x 


Bgure H-l. Method of Singularities Configuration 


The <r for supersonic flow past two oscillating blades 
in Zone II, (Figure H-1), can be obtained by superposition of 
Single blade solutions. For this purpose the disturbed flow 
field in Zone II is written as the linear superposition of the 
Single blade solutions for blades 1 and 2. However, in order 


to satisfy the flow tangency condition on the upper surface 


2) 5) 





of blade 2, a third single blade solution with unknown normal 


velocity must be added. Hence, one can write 


gp a EG aes (=>) 


meee! .xX-yvM--1 ae Gay 





Gia) 
. km2 
meee x) ~yy/M--1 ~lyery (*178) Ss a 
d5 Va =" vitije M 1 J 9 laz=z (xX, EF) (M l)y, ]daé 
4) 
. kM2 
a ee Xp-yy YMZ=1 — 1] ——— (x, -& ) kM a. ee 
orp a =" 9 Vv, (&)e M2eT 1 tS renee Oh E) (M l)y, de 


Vise) 

where yo denotes the unknown normal velocity. Equation (H-3) 
represents the disturbance due to the first blade, equation 
(H-4) the disturbance due to the second blade alone, and 
equation (H-5) the interference disturbance. Note that the 
blade 2 and interference potentials are expressed in terms of 
an Xy7 Yq coordinate system with origin at the leading edge 
Saeolade 2. 

Satisfying the boundary condition on the upper surface of 


blade 2 gives an equation for the unknown velocity, Vor tee; 





. kM? 
va (xryy= ee 2 PVM TD (ey ZT HE) 
WMZ-] “* 0 
Mk 
: Teac mea (Mo—l)y }dé 
(He=6)) 
at y=d. 
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Now setting 


YAS) Teale IVS ae lh 18), Osk= 1) 


and assuming the blades to pitch about their leading edges, 
eC. , 
Z(—E) = € ets)), 


Pmemetchnen performing the differentiation, one obtains 


eke 
Me txpy)= —— oY aks (6) 42° (Eye M1 8) 
g /Mz7=1 ) 0 
Toe he ee KM / (x=) 2- (Mey) de 


S/eeel| 
Cae a= iy 


. kM2 
- [ikZ (x-y /M7-1)+2~° (x-yYM7-1) ] ‘ itty, 
at y=da {H-9) 
Since, in general 


Cp = ikd + >. (H-10) 


the total pressure coefficient can be written as 


Cy = C51 a Ch5 = —- (H-11) 


where Col is obtained from the single blade solution for 


mrade | at y = d, 





, a - . a KM 
Beye (ta tc) 4827 (ee N21 FE) 
/iz=1 0 


Bo, (5% JETTA a6 -p 9 FO “Toa (ey +ha- Ce) ] 
km? 
e MZ-T 8). gM tre) = 1) a7 as 


Coe ee KM eZ (Mz= Taz 3 
+ ge °M?-T * JQ laz=y x M Te! Ib (H a2) 


sa) 





C59 from the single blade solution for blade 2 at a 0 


beetle a) 
Cp2 x ———(f mal [- 2(E) +227 (E) Je Mea 
—— 0 





; kM a eae ele wo eee ENC 
Moen (16) ce £48 fous (eb thao ey yo ao 





Pee a. -sepyar-+ 2 -i st Se ee (H-13) 
meet CU Ce Omer 1. , 
and 
site = ik dtp + orp (H-14) 
where 
l SoMa ay 
a a Lazo7 > mee 6 
eip — M ) olaeay Oxy - y) ]dy 
(H-15) 
and 
. kmZ 
a = IL Os _ kM? ~Ayzay (x oe) 
°rpx = — ie Wa 1az-1 e M il, 
kM 7 5 US 5 ea) kM 
: = = Son VICES = : 
Tq lorry (xy -y) Idy fe V Yen aaa al veal 
ee, (H-16) 
where . ; 
KM 
Oey = pe) ia (e)ee' (E) le tMz=1 (* 8) 
4 YM7=-1 0 
, RM ys OM gepe= terra ya 
¥(x-&) 2=(M“-1) d2 
a kM? a (H-17) 
—-[ikZ(x-dVYM2-1) + 27 (x-dYM7-1) Je ~“M¢-1 


Note that the substitution of Vv, (x) into equations .(H-15) and 


(H-16) for oop and ® I px requires one to relate x and y by 
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xX = y + d tan Bs (H-18) 
(see Figure H-1l) 
These integrals were computed using the Trapexoidal Rule 
and IBM subroutine BESJ. A listing of the program follows 


in Appendix I. 
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APPENDIX I 


PROGRAM LISTINGS 


ee OO) OO OO OO0O00 OOO0O00O DODO OOCODOCDOOO OOO OOO O0O0O0000°0 
AA QI TI FLD P= ODO SATAY SUV OE ODO HAA MPLIADM DDNOAN OT AO OHO ANMSFUVOM 
OOOO OO OOO 6 wed ed ed ht HOI OI OI NI OI OIA OI ON) O09 OO HO OO OF OT SPSS ST SS 
2 OO OO OO COO OOC OO OOO OO OCOOO OOO OO OD COO OOOOCO OO 000000 
eee OO OO O OOO OO QOODOOVODOOOO OOOO OO COO OO OOOC OO O0O00O00 
Pe OO OOO OOO CO OQCOD COO CVOOSV OO CDWDOVOOVOOO OOO OVD O00000O0 
OO OOOO OCOODOSC DODO COOD OO OCOSO CODCOD OOO OOOOOC OO 0O00000 
PI) DOOD ODD DOODO OO OOO OCOCOCO OD OOOOOOO OOO O0000 000000 
o~ 


~ ~ [ad ~< 
NI LJ oO © 
I ale rd Ne a 
iL aw! = ee — 
ee (9) WJ 56 NOS 
— 2a OW »= © 
~ Ce Cy UN 
at pf LW ~ On 
a o” bd mt 
oa ec _— i ba om 
= ae 6x) * Y) J 
Lu —a = wa bt J 
ae we We a ww 
VY) am Qw QA. EB a 
oe ~ re Y) 
raed | ~ & ore x ) 
<< O--Y 4 ae <I _) 
WO oor — OY <f 
—_ TOW oe ~ IO 
‘a we wb N Ir 
mao bet meer 1) mo C3 Ce 
Co LCAOD- 28 Ge ES 
— ro Ow) es Oo LU Higa fas 
b— COdwn — Cre OO 
WY) > ~<I- ~ J a 
mi) mS } — OQ > ra 
(oa N <io ol iF a Le 
uwjcD ~ & &L) 3B BE ~ 
he 2 Onmkbh eS) a b— <9 Co 
Wr aC) = a LL Lj 
<}-— TOBZN ~ ~ = aw 
cae Yur & a Oat Zz 
aT: CO ees C) oy es) uj oO 
> hom AW «& “ em co (ag 
=> ~ ©) ee) UL 1 bt Pag 
UL mma OO bk Yr> ~ 
Oat NOk = ~ <I <I 4 ‘) 
~O JILL ge uj <t uw é 
AK Ofrw « Me Pal Ce aN) 
On Owe m OC mf ue ~ 3 
r<t Fa «WW ~ 2 JV) = mh = ~ 
-—- OW ON — << ~- < uy SJ res = 
a) rm J LL 8 eo — SS eet we. SCD ys CD a8, 
= = oe Ge <~ ™ ot Oem2woaZz ae 
O AW ef ~ << Os“ Yew eTaitiv<cd call 
2) SeYdTe Zr SMM ~-oOn a J AOE KMOD w= 
Uj LL MW eXayD wr =) CS iN — © i<ai—st aS ON 
NI meee © wT © ~ Lo <p ty OO NEY me YY) Zz 
—C) MOTD JOY Oa “NO Jat H$NOALOYCaAS 
Wea A eOO*r® 2H A WW NU ect Cu ie HH OOONOF+ 
ae LUO Se TAO o ~- Nee O “O] << TeH# Sz tHMNOF 
wo ee ewctheVUYS— WL at NZ f re fk Wt ext kM YY J 
ZY cH SSOOO m4 OO Om -OZ— OD YW weqeNZSOAINOX WM 
mY OY eect azitt afte YD Ie 2S OFAZwwZszraze az 
Ul | Tei XeeDTeoa wWw Ogres Yaa -—- =O CO ewer WU a HY 
Oo 200 NSU tO SF aMe We DAN ODO hy) eee 
a) OO *eL *& ee <I arene O LL O Q. ea-atl SY MANxXKO ew Ll = CQ] a4 a 
Y — RBMOoOoCY «KX 2H WN ~Ow~~ 2 ew W OOH Hl i ux< = it WoW a 
Mzwaoww ww be bY) hm UNE oe CAV TT il acd col a0 fn es | co Ce Em age wg | 
ZOmeaPYOO 1 & —I<l wilw —I FD Onrnwsoitidcdeawhys OF WOOWFE 
WS ey eu SF ASST Orere OO OF OD HVRTtTOZOROR HO rb >a 
Y SSZz KMESRAIDS es Aro ft aK = be ee LEH eC EET 
a Cx maim my i UOOWwOow Ww WO OP ata anWnowe¢r WwW Ww LO 
uw OwWr mY A Cuuwe § SF YK Ch YO KR OTEK DOAOQTZErFe YH ws Dey 
at a ONTOS) aN 
a 
i} 
On ©) 
ond AI 
OO UW > 1 CS WOW STS TS) WOOO 
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